Anisotropic Dirac fermions in a Bi square net of SrMnBi2 
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We report the highly anisotropic Dirac fermions in a Bi square net of SrMnBi2, based on a first 
principle calculation, angle resolved photoemission spectroscopy, and quantum oscillations for high- 
quality single crystals. We found that the Dirac dispersion is generally induced in the (SrBi)^ layer 
containing a double-sized Bi square net. In contrast to the commonly observed isotropic Dirac cone, 
the Dirac cone in SrMnBi2 is highly anisotropic with a large momentum-dependent disparity of 
Fermi velocities of ~ 8. These findings demonstrate that a Bi square net, a common building block 
of various layered pnictides, provide a new platform that hosts highly anisotropic Dirac fermions. 

PACS numbers: 71.20.-b, 71.20.Ps, 71.18.+y, 72.20.My 



Unlike commonly observed quadratic dispersions, a 
linear energy dispersion, similar to the spectrum of rela- 
tivistic Dirac particles, is found in the so-called Dirac 
materials such as graphene and topological insulators 
(TIs). Such a linear dispersion is formed when two bands 
cross each other without hybridization due to the oppo- 
site spins or pseudo-spins. Linearly extended bands in 
Dirac cones, along with the (pseudo)spins, determine the 
peculiar properties of Dirac materials. For example, back 
scattering is strongly suppressed leading to a high elec- 
tron mobility, and anomalous half-integer quantum Hall 
effects are observed with a nonzero Berry's phase. [H, Hj. 

While the Dirac cones in graphene and TIs are in 
general isotropic ones, different types of Dirac cones 
have been proposed recently. Highly anisotropic Dirac 
cone was theoretically predicted in graphene under exter- 
nal periodic potentials [3] or mechanical stress [4]. Tilted 
Dirac cones with angle-dependent Fermi velocity have 
been suggested in a layered organic conductor a-(BEDT- 
TTF)2l3 at high pressures[5]. As a hybrid case, a semi- 
Dirac cone with a quadratic dispersion in one direc- 
tion and linear in the other may appear in hypotheti- 
cal graphene with a particular set of hopping parame- 
ters between the nearest neighbors [6]. Presence of such 
a hybrid Dirac cone has been proposed in quantum con- 
fined V02/Ti02 nanostructures.[7] While the realization 
of these Dirac fermions may lead to a discovery of novel 
electronic states e.g. di new- type of quantum Hall state [sl, 
there has not been a direct experimental observation of 
above mentioned Dirac Fermions. 

In this Letter, we report on a new bulk Dirac material 
SrMnBi2 that hosts highly anisotropic Dirac fermions in 
its Bi square net. Based on results from a first principles 
calculation, quantum oscillations and angle resolved pho- 
toelectron spectroscopy (ARPES) on high-quality single 
crystals, we demonstrate that there is a Dirac dispersion 



in the electronic structure of the double-sized Bi square 
net. The Dirac electrons have pseudo-spins stemming 
from wavefunction amplitudes at two different sublat- 
tice, similar to the case of graphene. In contrast to the 
graphene case, however, the Dirac cone in SrMnBi2 is 
highly anisotropic showing a significant momentum de- 
pendence in the Fermi velocity {vp) with a ratio between 
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FIG. 1: (color online) (a) SrMnBi2 crystal structure. Two 
building blocks, (b) MnBi layer and (c) Bi square net are 
shown. Note the double-sized unit cell in the Bi square net. 
(d) X-ray diffraction pattern from basal planes of a cleaved 
SrMnBi2 crystal, showing only the (00/) reflections. An opti- 
cal image of a single crystal with a scale of 250 /xm is shown 
in the inset, (e) Magnetic susceptibilities x{T) under H = 1 
T applied parallel and perpendicular to c-axis in both zero- 
field-cooled (ZFC, open) and field-cooled (FC, solid) runs. A 
kink in FC x{T) at Tn ^ 290 K is clearly seen in the inset, (f) 
In-plane resistivity (p(T)). It has a quadratic T-dependence 
at low temperatures as shown in the inset. 



the maximum and minimum vp^s of ~ 8. These find- 
ings suggest that the Bi square net, a common building 
blocks of various layered compounds [9], can provide a 
new platform for highly anisotropic Dirac fermions. 

Single crystals of SrMnBi2 were grown by melting sto- 
ichiometric mixtures of Sr (99.99%), Mn (99.9%) and Bi 
(99.999%) chunks in a sealed quartz ampoule. The am- 
poule was heat-treated at 1150 °C, followed by a slow 
cooling at a rate of 2 °C/h. X-ray diffraction and energy 
dispersive spectroscopy confirm high crystallinity[10] (see 
Fig. 1(d)) and right stoichiometry. Magnetotransport 
properties were measured in a conventional 6-probe con- 
figuration in high magnetic fields up to static 33 T at 
the National High Magnetic Field Laboratory(NHMFL) 
and up to pulsed 63 T at the Dresden High Magnetic 
Field Laboratory (HLD). ARPES data were taken with 
24 eV photons at the beamline 7U of UVSOR. Samples 
were cleaved in situ in an ultra high vacuum better than 
7 X 10~^^ Torr. First principles calculations were done 
using the full-potential linearized augmented plane wave 
method[l3] implemented in WIEN2k package[14]. The 
generalized gradient approximation (GGA) was utilized 
for the exchange correlation potential^]. 

We first look at the basic electronic properties of 
SrMnBi2. SrMnBi2 consists of a MnBi layer with edge- 
sharing MnBi4 tetrahedrons and a 2-dimensional (2D) Bi 
square net stacked with Sr atoms as shown in Fig. 1(a). 
Due to the low electronegativity of Sr, inserted Sr layers 
donate charges to the adjacent layers and also electroni- 
cally separate the MnBi layers and the Bi square net. In 
the [MnBi]- layer, Mn^^ has a half-filled 3d sheh in 3d^ 
configuration. The strong Hund couphng of Mn^+ leads 
to a magnetic ground state in MnPn {Pn = pnictogen) 
layers. [16j Magnetic susceptibility x(^) f^ct exhibits a 
kink at Tn ^ 290 K (Fig. 1(e)), indicating an antiferro- 
magnetic (AFM) transition. With the AFM ordering, the 
charge conduction in MnBi layers is highly suppressed. 
On the other hand, the covalent nature of Bi 6p bonds in 
the 2D square net makes the plane conducting, thus gov- 
erning the transport properties. Resistivity p{T) in Fig. 
1(f) indeed shows a highly metallic behavior, that is, a 
quadratic temperature dependence of p{T) = po + AT^. 
The parameter A, which is inversely proportional to the 
Fermi temperature, is found to be 19(1) nl^cmK"^, com- 



parable to that of pure Bi[ll|. This suggests that light 
carriers are responsible for the metallic conduction. 

The electronic structure of SrMnBi2 from a first prin- 
ciples calculation clearly reveals the origin of the light 
carriers. First of all, from the comparison of the to- 
tal energies for several magnetic states, we find that a 
Neel-ordered state in the MnBi layer is the most sta- 
ble state [ig]. With a Neel-ordering in the MnBi layer, 
we calculate the band structures of SrMnBi2 with the 
spin-orbit coupling(SOC) taken into account for the Bi 
bands (Fig. 2(a)). Due to the large spin polarization 
of Mn 3d electrons, Mn bands are placed away from the 
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FIG. 2: (color online) (a) Band structures of SrMnBi2, (b) 
total (gray line) and local Mn (red solid) DOS are shown in 
the upper panel. The lower panel shows local DOSs at Bi (red 
solid) and Sr (blue dashed) sites in the SrBi layer, (c) Band 
structure of the isolated (SrBi)+ layer where the line thick- 
ness represents the Bi 6px,y orbital character, (d) Anisotropic 
energy surfaces around the Dirac point ko = (0.208, 0.208) for 
the (SrBi)^ layer, (e) Dispersions near the Dirac point par- 
allel or perpendicular to the F-M symmetry line. SOC in the 
Bi bands is turned on for the results in (a) and (b), while 
SOC is not considered for the results in (c), (d) and (e). 



Fermi level {Ef) as seen in the density of states (DOS) 
in Fig. 2(b). The states near Ep are dominated by the 
Bi states in the square net where the Dirac-like energy 
dispersion is clearly seen at /co = (0.208, 0.208). A small 
gap near the Dirac point appears only when the SOC is 
taken into account, suggesting that this gap is not due to 
the orbital hybridization but from the SOC. Therefore, 
the light carriers are, in fact, Dirac fermions residing in 
the 2D Bi square net. 

In order to investigate the nature of the Dirac fermions 
in the Bi square net, we calculate the band structure 
of a single (SrBi)+ layer composed of a Bi square net 
stacked with Sr atoms above and below, and plot it in 
Fig. 2 (c). For the sake of clarity, we did not include 
the SOC in this case. We note that Dirac cones exist at 
several /c-points but only one at A:o=(0.208, 0.208) retains 
the Dirac dispersion in the band structure of SrMnBi2. 
Therefore, we focus on the Dirac cone at ko. Near the 
Ef^ states have mainly Bi 6px,y character with a small 
contribution from Sr Ad orbitals due to hybridization. Sr 
atoms below and above the square net causes unit cell 
doubling, resulting in two Bi atoms per unit cell. This 
leads to a folding of the dispersive Bi 6p bands. The 
two Bi 6px,y bands from two Bi sites cross each other 
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at a single point. There, the pseudo-spin of the Dirac 
cone in SrMnBi2 can be defined as the ampHtude of the 
wave- function at different Bi sites, similar to graphene. 

In spite of the similarity to the graphene, there are dif- 
ferences. First of all, the Dirac cone SrMnBi2 is highly 
anisotropic as seen in the Dirac dispersion of the (SrBi) + 
layer in Fig. 2(d). Along the T-M symmetry line, the 
Fermi velocity is = 1.51x10^ m/s which is compara- 
ble to that of graphene. In contrast, the dispersion along 
the cut perpendicular to T-M line is much weaker with 
Vp of ~ 1.91x10^ m/s, a factor of ~ 8 reduction com- 
pared to v^jp. The anisotropic nature of the Dirac cone 
is firmly retained in the full band structure of SrMnBi2. 
Existence of such anisotropy stems from the fact that the 
dispersion along the T-M line is determined by overlap 
between the neighboring Bi atoms in a square net while 
that along the line perpendicular to the T-M line is due 
to the hybridization between Sr dxz,yz orbitals and Bi 
Px^y orbitals. Therefore, the Dirac dispersion in a double- 
sized Bi square net is always highly anisotropic. Another 
difference is in the SOC gap. Unlike the negligible SOC 
gap in graphene[l8|, there is a sizable SOC gap of ^ 40 
meV at the Dirac point in SrMnBi2. This implies that, by 
placing the Fermi level inside the SOC gap, one can make 
an insulating Dirac system in which a quantum spin Hall 
effect is expected to appear [19]. These differences clearly 
demonstrate that the Dirac cone in the Bi square net of 
SrMnBi2 is very distinct from that in graphene. 

In order to experimentally investigate the elec- 
tronic structure, we performed ARPES experiments on 
SrMnBi2 single crystals. Figure 3(a) shows a FS map 
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FIG. 3: (color online) (a) ARPES intensity map at Ef- (b) 
Calculated FS. (c) Density map and (d) MDCs of the ARPES 
data along the T-M line (cut 1). (e) and (f) ARPES data 
along the cut perpendicular to the T-M line (cut 2). The cut 
directions are shown in panel (a). 



which is obtained by integrating the ARPES intensity 
within 25 meV energy window about the Ep. There are 
two distinct regions with high photoemission intensities: 
Large circular FS centered at the F point and a needle- 
like one between the F and M points. To compare the 
data with the calculation, we plot calculated FS's in Fig. 
3(b). There are three small FS's: Needle-shape FS's on 
the T-M line and near the X point, and fiower-shaped 
hole pockets near the Z point. In comparison with the 
calculated FS's in Fig. 3(b), the large FS near the F 
point should come from the fiower-shaped pockets near 
the Z point as the momentum resolution along the kz di- 
rection is relatively poor due to the finite escape depth of 
the photoelectrons. The needle-like FS on the T-M line, 
on which we focus, is consistent with the calculated FS 
shown in Fig. 3(b). On the other hand, the calculated 
FS's near the X point, originated from a band bottom, 
are not observed. The discrepancy is most likely due to 
the underestimation of the band gap near X-point in the 
GGA exchange correlation functional. 

Now we focus on the needle-shaped FS on the T-M line 
where a Dirac dispersion is expected. Figs. 3(c) and 3(e) 
are ARPES intensity plots along the cuts parallel and 
perpendicular to the T-M line as indicated in Fig. 3(a). 
From the raw data and the corresponding momentum 
distribution curves (MDCs) depicted in Figs. 3(d) and 
3(f), energy dispersion with a strong anisotropy is iden- 
tified. For the cut 1, we mark the dispersion determined 
from the MDCs in figure 3(d). The experimental band 
is almost linear at higher binding energies but appear to 
show some curvature near the Ef. Such behavior is con- 
sistent with the calculated band with a shift of chemical 
potential ~ 0.08 eV. in Fig. 3(c) which shows an SOC 
driven gap. The Fermi velocity I'i? is ~ 2 x 10^ m/s for 
the cut 1. In a stark contrast, the band for the cut 2 in 
Fig. 3(e) is extremely dispersive with a much larger vp 
(> 1 X 10^ m/s). In spite of the fact that presence of the 
SOC does not allow a Dirac point, the overall dispersion 
shows a Dirac nature with a strong anisotropy. 

As we have identified an anisotropic linear dispersion 
at a high energy scale, it will be important to investi- 
gate its chiral nature due to pseudo-spins. For that pur- 
pose, we have performed magnetotransprot experiments 
in high magnetic fields. Fig. 4(a) shows the in-plane 
magnetoresistivity (Ap^x) and Hall resistivity (pxy) as 
a function of out-of-plane magnetic field. Clear quan- 
tum oscillations in both pxx{H) and pxy{H) are observed. 
The background pxx{H) exhibits a linear H dependence 
without any signature of saturation up to 63 T. Such a 
non-saturating linear magnetoresistance has been often 
found when the conditions, no <C (eB. jhcf"!^ and m* <C 
eRhjcT^ are satisfied [l2|. The crossover magnetic field 
between the regions of Apxx ^ and Apxx ^ H is only 
~ 1 T. This implies that the charge conduction is deter- 
mined by low-density carriers with an extremely small 
effective mass, which is consistent with the Shubnikov- 
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FIG. 4: (color online) (a) Magnetic field dependence of 
Apxx{H)/ po and pxy of SrMnBi2 (b) Hall SdH oscillations 
A(Jxy = (Jxy - {(^xy) as a function of 1/H at various temper- 
atures, where (Txy — pxy/{pxx~^Pxy)- (c) Cyclotrou mass and 
dingle plots for Aaxy. (d) SdH fan diagram plotting the mea- 
sured 1/Bn with the filling factor n which is estimated from 
the Aaxy versus 1/Bn plot shown in the inset. 



de-Hass(SdH) oscillation results in Fig. 4(a). The oscil- 
lating component shows a periodic behavior in 1/5 with 
a single frequency of F = 152(5) T. Using the Onsager 
relation F = {^o/27r'^)Ak^ where is the flux quantum 
and Ak is the cross-sectional area of the Fermi surface 
normal to the magnetic fields, Ak is found to be 1.45(5) 
nm~^ which corresponds to the ^ 0.7% of the total area 
of the Brillouin zone. The H or T dependence of the os- 
cillating amplitude (Fig. 4(b)) gives a cyclotron mass of 
rric = 0.29(2) rrie {rrie = free electron mass) and scatter- 
ing time of r = 3.5(5) x 10~^^ sec~^. The corresponding 
mobility /i = er/rric is ~ 250 cm^/Vs. A small FS, small 
effective mass, and relatively large mobility are consistent 
with the presence of Dirac fermions. 

The key evidence for the existence of Dirac fermions 
is the non-zero Berry's phase associated with their cy- 
clotron motion. According to a semi-classical magneto- 
oscillation description, the oscillation part of pxx follows 
^Pxx ^ cos[27r(F/5 + 1/2 + /3] where /3 is the Berry's 
phase (0 < /3 < 1). Experimentally, the Berry's phase 
can be obtained from the analysis of the SdH fan dia- 
gram plotting 1/Bn as a function of the Landau index 
n which can be determined from the minima of pxx as 
shown in the inset of Fig. 4 (d)[20j. In the SdH fan dia- 
gram, the intercept of the linear fit yields Berry's phase 
which is expected to be 1/2 for Dirac fermions in e.g. 
monolayer graphene. As shown in Fig. 4(d), the linear 



fit does not go through the origin but rather intercepts 
the the abscissa at n = -0.40(9). The resulting non-zero 
Berry's phase reveal that the low-density and high mobile 
carriers in SrMnBi2 are indeed Dirac fermions. 

Based on our findings, we conclude that the Bi square 
net in SrMnBi2 hosts highly anisotropic Dirac fermions. 
This implies that SrMnBi2-type compounds containing 
a Bi square net can be a new platform for anisotropic 
Dirac fermions with intriguing properties. For example, 
in AeTPn2 {Ae = alkaline earths, T = transition met- 
als and Pn = pnictogens), the anisotropy of the Dirac 
cone can be tuned by introducing other alkaline earths, 
and the SOC gap can also be tuned by replacing Bi with 
other pnictogens with lower atomic numbers. In addi- 
tion, since the magnetic ordering of the MnBi layers ad- 
jacent to the Bi square net is intimately related to isola- 
tion of the Dirac cone from other bands, one could con- 
trol the coupling between AFM and Dirac fermions by 
replacing Mn atoms with other transition metals hence 
modifying the magnetic order. Furthermore, supercon- 
ductivity can also be incorporated into materials with 
Dirac fermions as found in a recently-discovered super- 
conductor CeNia;Bi2[21] and a hypothetical supercon- 
ductor BaFePn2[22]. Therefore, our findings on Dirac 
fermions in a layered pnictide, SrMnBi2, containing a 
square net provide new perspectives not only for control- 
ling the anisotropy of the Dirac cone but also for coupling 
to magnetic or superconducting orders. 
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